Diffusion-limited unbinding of small peptides from PDZ domains 
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PDZ domains are typical examples of binding motifs mediating the formation of protein-protein 
assemblies in many different cells. A quantitative characterization of the mechanisms intertwining 
structure, chemistry and dynamics with the PDZ function represent a challenge in molecular biol- 
ogy. Here we investigated the influence of native state topology on the thermodynamics and the 
dissociation kinetics for a complex PDZ-peptide via Molecular Dynamics simulations based on a 
coarse-grained description of PDZ domains. Our native-centric approach neglects chemical details 
but incorporates the basic structural information to reproduce the protein functional dynamics as it 
couples to the binding. We found that at physiological temperatures the unbinding of a peptide from 
the PDZ domain becomes increasingly diffusive rather than thermally activated, as a consequence 
of the significant reduction of the free energy barrier with temperature. In turn, this results in a 
significant slowing down of the process of two orders of magnitude with respect to the conventional 
Arrhenius extrapolation from low temperature calculations. Finally, a detailed analysis of a typical 
unbinding event based on the rupture times of single peptide-PDZ contacts allows to shed further 
light on the dissociation mechanism and to elaborate a coherent picture of the relation between 
function and dynamics in PDZ domains. 



I. INTRODUCTION 

The role of PDZ domains in the organization of protein 
complexes at the plasma membrane has been increasingly 
recognized in the last decade^. Proteins containing sev- 
eral PDZ domains (up to 13 in the MUPP1 protein 2 .) 
act as scaffolds that cluster together different transmem- 
brane, membrane associated and periplasmic proteins in- 
volved, among other functions, in signaling pathways—^ 
and ion permeability 5 . The participation of PDZ do- 
mains in the organization of supramolecular complexes 
in skeletal muscle cells has also been documented 6 '?. 

PDZ domains associate with other proteins by binding 
their carboxyl-terminal aminoacids^, as highlighted by 
the structure of several PDZ domains, such as the third 
PDZ domain of PSD95 (Postsynaptic density-95/disks 
large/zonula occludens-1) 10 - and the PDZ domain of 
ZASP (Z-band alternatively spliced PDZ-motif) li, al- 
though internal structures, such as /3-hairpins that mimic 
C-terminal geometries can also be recognized, as in the 
complex between nNOS (neuronal nitric oxide synthase) 
and syntrophini^. 

Recently, due to their central role as key mediators 
of protein-protein interactions in mammalian cells, PDZ 
domains have been the object of intense study, with the 
aim of designing small molecules capable of acting as 
modulators or inhibitors of the PDZ binding activity in 
a controlled fashion. Efforts have focused in the design 
of both non-peptid o 13 i 14 ' 15 and peptid o 16 i 17 ' 18 ligands, 
with the ambition to develop molecular probes to study 
the biophysical and biochemical properties of PDZ do- 
mains and to devise new small molecule-based thera- 
peutic strategies. Hence, the great importance of un- 
derstanding the principles of peptide-PDZ interactions. 

The geometry and chemistry of binding to PDZ do- 
mains involve the fit of the last 4 to 5 carboxyl-terminal 
aminoacids into a groove between a cv-helix and a /3- 



strand on the PDZ surface (FigOJ, with the last C- 
tcrminal residue almost invariably hydrophobic. The 
specificity of each domain is conferred by a few (2 to 
3) PDZ surface aminoacids that make contacts with the 
residues in positions -1 to -4 relative to the C-terminal 
in the target proteiiA The surprising simplicity of this 
binding scheme possibly explains why PDZ domains are 
one of the most widespread binding modules yet identi- 
fied, since just a few incremental, concerted mutations in- 
volving surface aminoacids (hence unlikely to change the 
overall protein stability) can tune the affinity of PDZ do- 
mains for different targets. On the other hand, a binding 
architecture that relies on just a few optimized contacts 
comes at the price of losing strict specificity. Indeed, re- 
cent experiments on 26 mouse PDZ domains and domain 
clusters have confirmed that each PDZ domain can bind 
to several peptides^, and that each peptide, in turn, can 
bind to several PDZ domains. 

In a recent work 20 , we have explored via normal mode 
analysis (NMA) the mechanical aspects of binding of 
a peptide to a PDZ domain. In line with the results 
of NMA analysis of several protein o 21 ' 22 ? 23 , the picture 
that emerged was that a limited number of low- frequency 
modes suffice to reconstruct the observed conformational 
change from the apo to the complexed structure of the 
PDZ. The long-range spatial correlations that character- 
ize these modes correspond to a concerted breathing mo- 
tion of the binding cleft, thus suggesting that functional 
dynamics is deeply rooted in the native architecture it- 
self. However, despite their success in providing a qual- 
itative picture of the coupling between thermal fluctua- 
tions of the structure and the binding deformation, the 
predictive power of normal modes is still limited by the 
harmonic approximation. Hence, in order to investigate 
the progressive detachment of the peptide from the PDZ 
structure, one must resort to more general computational 
studies. 
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Molecular dynamics simulations have played a crucial 
role in understanding the basis of the concomitant se- 
lectivity and promiscuity of the PDZ binding dynam- 
ic o 24 i 25 . Basdevant et al£^, in particular, have used all- 
atom Molecular Dynamics (MD) with a realistic force- 
field to correctly reproduce the experimental ranking, 
but not the precise values, of binding free-energies. The 
chemical details of the target peptide and of the binding 
groove have been taken into account, and hydrophobic 
effects have been shown to be key determinants of the 
promiscuity, with other interactions conferring stronger 
or weaker selectivity. 

The free energy difference AG^ between the unbound 
and bound states is customarily related to the dissocia- 
tion constant Kd by Kd = exp(— AG^//csT) where T is 
the absolute temperature and ks is the Boltzmann con- 
stant. However, the constant Kd is also related to the 
dynamics of the system by the equation Kd = k s/k on , 
where k g and fc on are the unbinding and binding rates 
respectively. At a first approximation, k s and fc on de- 
pend on the free energy difference between the bound and 
unbound states and the barrier between them. The dis- 
sociation constant Kd, therefore, conceals the kinetics of 
the system, the same value being compatible with both 
very fast and very slow rates, as long as their ratio does 
not change. In the case of PDZ domains, dissociation 
rates range experimentally 26,27 from a few s _1 to 10~ 3 
s _1 , thus implying that a correct description of the un- 
binding process would need to capture the PDZ dynamics 
over comparable timescales. Unfortunately, detailed all- 
atom MD simulations can at present cover a few tens 2 ^ 
and at most reach a hundred nanoseconds. Therefore, in 
order to observe a typical unbinding event, MD simula- 
tions should be from 6 to 10 orders of magnitude longer 
than currently possible. Hence, one must resort to sim- 
pler, coarse-grained models in order to study the kinetics 
of interaction between a PDZ domain and its target pep- 
tide. 

Neglecting the atomic detail clearly has both advan- 
tages and drawbacks. On the one hand, the reduction in 
the number of degrees of freedom and the simplified force- 
field of interaction allow to explore longer times scales, 
thus making a complete equilibrium description of the 
binding kinetics possible. On the other hand, the chem- 
ical specificities can be reincorporated in the model only 
at a qualitative level. 

In this paper, we employ the Go strategy 2 ^, a native- 
centric scheme, coarse-grained at the residue level to 
simulate the unbinding dynamics of a peptide from the 
third PDZ domain of PSD-95 (henceforth referred to as 
PDZ3). Such simplified description has already proved 
successful in the characterization of the role of the native 
state topology and of its dynamics in protein recognition 
and binding mechanism s 29 ! 30 . Our simplified simulation 
scheme allows us to draw a clear, though approximated, 
picture of the dissociation kinetics. 

The escape of the peptide from the binding groove on 
the PDZ surface takes place over a free-energy landscape 




FIG. 1: Cartoon representation of the crystal structure of the 
PDZ3 domain in complex with its target peptide (blue). PDB 
code 1BFE. 

that strongly depends on the temperature. At physiolog- 
ically relevant temperatures, less than 10% lower than 
the unfolding temperature (measured to be about 320K 
for the second PDZ domain of PTP-BL 3 ^), the free en- 
ergy difference between the bound state and the barrier 
turns out to be rather small (0.5 kcal/mol) because of the 
nearly complete compensation of the enthalpic and en- 
tropic components. As a consequence, unbinding at phys- 
iological temperatures is eminently a diffusive, rather 
than thermally activated, process and the unbinding rate 
is orders of magnitude smaller than an Arrhenius-like ex- 
trapolation from low temperatures would suggest. 

II. METHODS 

The PDZ3 domain (PDB codes 1BFE and 1BE9 with- 
out and with bound peptide respectively) as resolved by 
X-ray crystallography^ is 110 residue long but we trun- 
cated the chain from Arg309 to Ser393 because the final 
loop does not participate to the binding mechanism and 
its large flexibility may hinder some signals specifically 
involved in the binding. 

We adopt the model proposed by Clementi et al2&. 
where successive beads along the chain are connected by 
stiff harmonic springs, mimicking the peptide bond and 
maintaining the chain connectivity, 

Vbb(n,i+i) = ^Ke(ri^ l+1 - i? Vi+ i) 2 (1) 

with stiffness K = 1000/d 2 ,, where d = 3.8 A is the 
mean equilibrium distance of two consecutive residues 
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along the chain, and e = 0.84 kcal/mol (see Results) sets 
the energy scale. Here, and R^j indicate the distance 
between residues i and j in a generic conformation and 
in the native structure (1BFE and 1BE9), respectively. 

In line with native-centric schemes, non-bonded (nb) 
interactions between non consecutive a-carbons are mod- 
elled with Lennard- Jones 12-10 potentials if the atoms 
are in contact in the native state according to a given in- 
teraction cutoff R c and with purely repulsive interactions 
otherwise 
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The parameters of the non-bonded interactions are fixed 
as R c — 6.5 A and a = 4.5 A. The force field is completed 
by the angular interactions 
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where 6i is the bending angle identified by the three con- 
secutive C Q 's i — 1, i, i + 1 and <pi is the dihedral angle 
defined by the two adjacent planes formed by four con- 
secutive C Q 's at i — 2, i — 1, i, i + 1. The superscript 
identifies quantities referring to the native conformation. 
The force field parameters are proportional to the energy 
scale e so that kg = 20e, kf^F 1 = e, fc^ 3 ' = 0.5e and one 
time unit corresponds to about 3 ps, when considering 
an average aminoacid mass of 110 Da. 

We have performed fixed-temperature Molecular Dy- 
namics simulations within the isokinetic scheme^ 2 ., which 
provides a correct sampling of the configuration space. 
We have then applied the multiple histogram technique 33 
to estimate thermodynamic observables such as the in- 
ternal energy, the specific heat 
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and the structural similarity Q of a given conformation 
with the native structure. The latter parameter is defined 
as 



E l>] Q(Rc-Rj)Q(Rc 



J2i >3 ®(Rc - Ri. 



(5) 



0(u) being the unitary step function. Q represents the 
fraction of native contacts present in a given conforma- 
tion. The multiple histogram technique allows to con- 
struct the free energy profiles G(Q) = — ksT log P(Q) as 



functions of Q, which plays the role of a reaction coordi- 
nate. 

In order to obtain the free energy profiles of unbind- 
ing, we employed the umbrella sampling technique^. In 
particular, we restrained the distance between the pep- 
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tide and the PDZ3 centers of mass p — \R p dz 
a given range of values (k = 1, 2, . . . , iV samp ) via the 
harmonic umbrella potential 

with K q — 0.84 kcal/mol/A 2 . The multiple histogram 
technique allows then a de-biasing at a given tempera- 
ture by matching all the histograms for the center-mass 
distance collected around the different sampling values. 
By varying the reference temperature we can obtain dif- 
ferent free energy curves. 

The statistics of the unbinding times have been col- 
lected through Langevin dynamics. The unbinding time 
is defined here as the time at which the bond distances 
of all the 13 contacts linking the peptide to the PDZ3 
exceed a given threshold for the first time. We fixed 
such thresholds at 1.5 times the corresponding values in 
the native conformation. However, small variations of 
this threshold did not result in major changes of the exit 
times statistics. 



III. RESULTS 

The energy scale e can be fixed by comparing the 
folding temperature X 1 / of the model, in units of e/ks, 
with experiments. For this purpose we studied the PDZ- 
domain specific heat C v as a function of temperature. 
The thermogram in Fig. [2] allows to identify the unfold- 
ing transition at T/ w 0.76. From available data 3 ^, the 
folding temperature of PDZ3 is close to 323 K (50° C), 
which allows to set the model energy scale to the value 
e = 0.84 kcal/mol. The free energy profile as a func- 
tion of the fraction of native contacts Q (Fig. [21 inset) is 
typical of clean two-state transitions, in agreement with 
experiments 3 ^, where folding was found to be affected at 
most by a high energy poorly populated intermediate. 

A further check of the chosen energy scale may be ob- 
tained from a comparison with the experimental stabil- 
ity AG = 6-6.5 kcal/mol of PDZ3 at 298 rOl. From 
our simulations we can estimate Go as the energy differ- 
ence between the minima in the native and denaturated 
basins. At T — 298 K we find AG = 5 kcal/mol, in 
good agreement with experiments. 

It is our aim to investigate the unbinding dynamics of 
a small peptide from the PDZ3 domain. Hence, both 
the unbinding and unfolding transitions have to be care- 
fully located in temperature. The free-energy profiles as 
a function of the distance between the centers of mass of 
the protein and of the peptide, chosen as reaction coordi- 
nate for the unbinding, are reported in Fig. [3] for several 



4 




200 300 400 50 

r(K) 

FIG. 2: Specific heat curve of Go-model PDZ domain fold- 
ing/unfolding. The inset shows the double well feature of the 
free energy profile (kcal/mol) as a function of the fraction of 
native contacts Q at the folding temperature. 

temperatures lower than that of unbinding. At low tem- 
peratures, the curves are characterized by a well defined 
minimum corresponding to the bound state and by a 
steep free energy barrier that has to be overcome in order 
for the complex to thermally dissociate. Remarkably, the 
barrier turns out to decrease rapidly at increasing tem- 
peratures, marking an almost complete compensation of 
the enthalpic barrier by the increase of conformational 
entropy. The behaviour of the dissociation free energy 
AG(T) = G ban (T) - G min (T) with the temperature is 
shown in the inset of Fig. [3] As expected, AG is a de- 
creasing function of T, vanishing at Tdiss = 306 K, which 
locates the spontaneous unbinding 17° C below the un- 
folding temperature. 

A. The unbinding dynamics 

The data shown in Fig. [3] contain the relevant informa- 
tion on the unbinding dynamics. In fact, it is possible to 
compute the dissociation rates by straight integration of 
the numerical free energy profiles according to Kramers' 
theorjs^, 

fc ° ff ( T ) = fX* G(x) f X G(x') ^ 

7 / e k B T dx e k B T dx' 
Jx Jo 

where xq and Xb are the abscissas of the minimum and 
of the barrier maximum, respectively. The temperature 
curve of the dissociation rates as calculated from for- 
mula (|6|) is plotted in Fig. [4] Our results clearly show that 
the simple Arrhenius description k Q g cx exp(— AG/fc^T) 
only holds in the low temperature region, where the free 
energy barrier is high enough to justify the familiar treat- 
ment of a thermally activated unbinding. At higher tem- 
peratures the dissociation process slows down consider- 
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FIG. 3: Free energy profile along the reaction coordinate at 
different temperatures below the unfolding of the complex 
PDZ3-Peptide. From top to bottom, the curves refer to tem- 
peratures from T = 85 K up to T = 255 K in steps AT = 21 
K. The inset shows the temperature dependence of the free 
energy barrier (kcal/mol). 

ably with respect to the extrapolated low-temperature 
prediction. 

This effect is due to the gradual switch from a barrier- 
limited, thermally activated process at low temperatures 
to a diffusion-limited process as the unbinding temper- 
ature is approached. Indeed, using ([6]) in the limit case 
where there is a perfectly flat free-energy landscape be- 
tween xq and Xb (and a reflecting wall in xq), hence no 
barrier, k s = ^{ksT f^) / {xb — xq) 2 , which is the typical 
rate for particles diffusing from xq, with diffusion con- 
stant ksT/j, and absorbed at Xb- The diffusion-limited 
k e depends only slightly on temperature when repre- 
sented on an Arrhenius plot (Fig. QJ. 

However, the straightforward application of Kramers' 
theory at high temperatures may not be legitimate due 
to the strong reduction of the free energy barriers, that 
makes the hypotheses of local equilibrium in the free en- 
ergy minimum and high barrier questionable. 

In order to check the validity of the Kramers' predic- 
tion at high temperatures, we carried out kinetics simu- 
lations of unbinding to directly estimate the dissociation 
rate at a temperature T. We observed the time evo- 
lution of the PDZ3 /peptide complex until a spontaneous 
unbinding event was recorded. By doing so, we were able 
to measure very accurately the lifetimes r of the molecu- 
lar complex, the temperature-dependent average (r) and 
distributions P{t). The inverse of (r) provides a mea- 
sure of the kinetic dissociation rate: k a s = 1/(t). The 
results of this procedure are shown in Figure 0] together 
with the rates estimated from Kramers' formula ©. The 
agreement between unbinding simulations and the ana- 
lytical estimates is remarkable at all temperatures where 
we were able to collect enough statistics for the compu- 
tation of the average inverse lifetime of the complex. In 
fact, dissociation events become exponentially rarer at 
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FIG. 4: Arrhenius plot comparing the dissociation rates esti- 
mated from the free energy profiles of Fig. [3] through Eq. © 
and the rates computed as inverse average times from di- 
rect unbinding simulations. Since pre-factors are unknown 
by construction, the latter have been rescaled by an appro- 
priate temperature-independent factor to match the values 
obtained from Kramers' theory. All rates are normalized by 
the kinetic rate at room temperature T a — 300 K. The thick 
dashed line is the result of an Arrhenius fit with the expres- 
sion fc fj = a exp(— A/ksT), which gives a free energy barrier 
height of A ~ 7 kcal/mol . The vertical dashed line marks 
the room temperature. 




FIG. 5: Histogram of the kinetic unbinding times at T — 
300 K (symbols). The solid line in the inset is a fit of the 
early times region with the first passage time distribution of 
a diffusive process P(t) = <zexp(— tq/t). 



low temperatures, thus requiring unrealistic simulation 
times in order to collect satisfactory statistics. These 
results provide a direct confirmation of the validity of 
Eq. ^ in the high temperature region. 

That the unbinding process is diffusion-limited at high 
temperatures is nicely confirmed by the analysis of the 
the histogram of dissociation times. In Fig. [5] we show 
one of such distributions computed at T = 300 K. At 
long times P(t) decreases exponentially, implying that 
the slow unbinding events obey the statistics of the wait- 
ing times of a Poisson-like process, as should be expected 
for a simple Arrhenius picture of thermal activation over 



a barrier. However, it can be clearly appreciated that 
fast events are under-represented with respect to what 
the Arrhenius picture would predict. In particular (see 
inset of Fig. [SJ, we find that the left tail of the distri- 
bution follows the law P(t) oc exp(— tq/t), which indeed 
characterizes the distribution of first passage times in a 
diffusive process. 

As a final methodological remark, we wish to com- 
ment on using the Langevin scheme for the unbinding 
simulations. It is known that the kinetics simulated with 
Langevin MD is strongly dependent on the friction co- 
efficient 7. We thus performed different runs at differ- 
ent values of 7 in order to check the dependence of the 
unbinding kinetics on the friction. The results of this 
analysis are reported in Fig. [51 where we plot the av- 
erage peptide dissociation time (t) at an intermediate 
temperature T = 234 K as a function of 7. We explored 
both the under-damped and the over-damped regimes, 
and compared the numerics with the asymptotic limits 
of Kramers' theory of barrier crossing^ 
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where luq and u>b are the frequencies at the bottom and 
saddle points of the free energy profile and 
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y/-2m[G{x) - G(x b )} dx , 



x\ being the left-hand side turning point, i.e. the solu- 
tion of the equation G(x\) — G{xb) with x\ < xq. Wc 
find that our data indeed converge toward Kramers' esti- 
mates in both the under-damped and the over-damped 
regimes, thus showing that the computation of k Q s as 
1/(t) can be meaningfully extrapolated to different val- 
ues of 7 at different temperatures. In particular, we 
can make use of the theoretical prediction for the over- 
damped regime, seen to be valid for values of 7 approxi- 
mately greater than 0.1, for testing whether the results of 
Fig. [4]still hold for other values of the damping parameter 
when using frequencies and free energy barriers computed 
numerically. We see from Fig. [7] that the essential fea- 
tures of the Arrhenius plot are preserved. Importantly, 
the room temperature value of the dissociation rate re- 
mains still largely overestimated if extrapolated from the 
low-temperature values even at larger damping 7. 

The scenario outlined above is consistent with other 
MD studies employing the Langevin scheme to investi- 
gate the kinetics of folding—. 



B. Unbinding dynamics of individual contacts 

A more detailed description of the unbinding process 
can be obtained by considering separately the kinetics of 
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FIG. 6: Average kinetic rate at the intermediate tempera- 
ture T = 234 K as a function of the friction 7 used in the 
Langevin simulation of the unbinding kinetics. The solid and 
dashed lines are the predictions of Kramers' theory ([7} calcu- 
lated with the frequencies and free energy barrier estimated 
from the corresponding simulated free energy profiles of Fig. [3] 
Note that no fitting is performed here. 
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FIG. 7: Arrhenius plot comparing the dissociation rates esti- 
mated for different values of 7 from the over-damped formula 
in eq. Q with the temperature-dependent parameters u>o, cut 
and AG estimated from the simulated free energy profiles. All 
rates are normalized by the kinetic rate at room temperature 
T a — 300 K for 7 = 0.05. The vertical dashed line marks the 
room temperature. 



the 13 contacts that keep the peptide bound to the PDZ 
in the native state. Each of them will be characterized 
by a different breaking time (r c ), that can be estimated 
as the average over many dissociation events in a typical 
unbinding simulation. The ranking between such times 
provides an average picture of the order at which contacts 
break down during the gradual detachment of the peptide 
from the PDZ domain. 

It is thus interesting to investigate whether there ex- 
ists a correlation between the unbinding event of each 
PDZ/peptide bond and the dynamics of the residues 
from the PDZ domain involved in the same link. A rela- 



tion between the binding dynamics and a specific spatial 
pattern characterizing a reduced subset of low-frequency 
normal modes has been recently suggested by De Los 
Rios et al£Q. In the present case, we find four modes 
among the eigenstates of the Hessian of our force field 
that match the pattern identified in Ref. [2(| namely the 
set S — {7,8,9,14}. We thus introduce an indicator 
measuring the spectral weight corresponding to that set 
of modes at each site 



Wj 



E E & 

fcgS a=x,y,z 



(8) 



where (g a is the a (a — {x, y, z}) component of the fc-th 
normal mode at site i. In figure [5] we plot (t c ) vs. Wi for 
the 13 binding contacts, where i indicates the residue of 
the PDZ domain participating to the contact c with the 
peptide. 

The data show to be clustered into three sets. The 
right outermost cluster contains the two contacts that 
have both a short lifetime and the larger spectral weight. 
In terms of the latter indicator, these contacts are char- 
acterized by the fourth and fifth largest values among all 
residues. This strongly suggests that functional normal 
modes are likely to contribute substantially to the loos- 
ening of the PDZ /peptide bond and its eventual rupture. 
The left bottom cluster involves the contacts with small 
(t c ) but whose local fluctuations are poorly captured by 
the functional modes. Taken together, these first two 
clusters include the residues forming the hydrophobic 
pocket: the second group comprises the residues at the 
bottom of the pocket whereas the first one includes two 
residues flanking the entry of the binding pocket. Finally, 
the last contacts to break up involve the less conserved 
residues that confer specificity to the different PDZ do- 
mains. 

Overall, these results outline the following picture of a 
typical unbinding event. The first contacts to loosen up 
are those involving residues from the hydrophobic bind- 
ing pocket, part of which are involved in the breathing 
pattern of the cleft as highlighted from normal mode 
analysis. This observation provides a nice example of 
the interpretation of a specific normal mode pattern as 
precursory feature of a much larger amplitude motion. 
The contacts that break up the last are those from the 
surface region of the binding pocket, the ones associ- 
ated with binding specificity and thus the less conserved 
ones. Since these residues are obviously more important 
as their chemical details are concerned, it is not surprising 
that a dynamical investigation based on a purely topo- 
logical model does not assign them a key role in the un- 
binding kinetics. 



C. Conclusions and Discussion 

The central result of our analysis of the unbinding of 
peptides from the PDZ3 domain is that the bound state 
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Spectral weight of PDZ residue 



FIG. 8: Average breaking times of individual contacts ver- 
sus the spectral weight on the subset of functional modes of 
the corresponding residue in the PDZ domain. The breaking 
times are normalized to the average unbinding time of the 
peptide as a whole. 

is stabilized over a physiologically relevant temperature 
interval. 

Countcrintuitively, this happens as a direct conse- 
quence of the strong temperature-dependent quenching 
of the free energy barrier in the proximity of the func- 
tional region of temperatures. As a matter of fact, such 
slow down of the dissociation reflects the gradual change 
in nature of the fastest events from thermally activated 
to diffusive. Indeed, as the free energy barrier gets sub- 
stantially reduced at increasing temperatures due to the 
entropy-enthalpy compensation typical in the proximity 
of the unfolding transition, Kramers' formula ([6]) pre- 
dicts a diffusion dominated unbinding process, and ki- 
netic simulations confirm such picture. Overall, diffusion 
slows down the process of two orders of magnitude with 
respect to the prediction of the Arrhcnius formula ex- 
trapolated from the low-temperature region. 

The spectral analysis of the breaking time of the con- 
tacts involved in the binding helps shedding further light 
on the unbinding process, that on average first involves 
the escape from the hydrophobic binding pocket and sub- 
sequently from the surface region where residues confer- 
ring specificity are usually located. From the time his- 
tory of all the contacts that keep the peptide bound to 
the domain until they rupture we obtained their aver- 



age lifetimes, that correlate with a picture emerged form 
NMA analysis. This suggests a non trivial coupling be- 
tween collective dynamics of the molecule and binding 
mechanism even at physiological temperatures. This con- 
firms that Normal Mode Analysis may provide a valuable 
method for detecting the protein chain fluctuations that 
represent the precursory events to unbinding. 

Our results rest on the reliable calculation of the free- 
energy profile of unbinding, which is a quantity defined 
at (quasi-)equilibrium, and on the precise determination 
of the escape-time histograms, which implies a complete 
description of the metastable state of the process, again 
a (quasi-) equilibrium quantity. This is consistent with 
the observation that experimental unbinding rates range 
from 1 to 10~ 3 s _1 , time-scales long enough for nanomet- 
ric dynamical molecules such as PDZ domains to fully ex- 
plore their phase-space and consequently to be amenable 
to a (quasi-)equilibrium description. 

Unfortunately, present all-atom MD with realistic 
force-fields can access at most few hundreds nanoseconds 
of the dynamics, and are thus likely to miss events rele- 
vant on the second or longer time-scales. Employing a na- 
tive centric backbone representation of the protein with 
simplified force field (Go-like model) on the one hand al- 
lowed us a thorough investigation of the equilibrium and 
stationary properties of the system, and on the other 
hand emphasized the role of the native state geometry in 
the binding mechanism of the PDZ3 domain. Given the 
overall high structural similarity of the PDZ family, our 
evidence for a diffusion-limited unbinding process likely 
applies to most if not all of them. 

Although necessary to access the relevant time-scales, 
neglecting the chemical details of the binding cleft and 
of the peptide comes at a price: it is impossible to pre- 
cisely rank the affinities of different peptides for the same 
domain and of different domains for the same peptide. 
At the same time, the promiscuity of PDZ binding im- 
plies that some degree of generalized stickiness, captured 
by the simplified force-field used here, is present. More- 
over, we believe that our results represent a reading frame 
for future unbinding simulations that will become possi- 
ble when algorithms and computational resources will be 
powerful enough to access the experimental time-scales 
of the process while taking into account the full chemical 
detail. 

At the same time, we stress that our prediction of diffu- 
sive vs. thermally activated unbinding is surely amenable 
to experimental verification. 



1 Fanning, A. S., Anderson, J. M., J. Clin. Invest. 1999, 103, 
767-772. 

2 Ullmer, C, Schmuck, K., Figge, A., Labbert, H. FEES 
Lett. 1998, 424, 63-68. 

3 Pawson, T.; Scott, J. D. Science 1997, 278, 2075-2080. 



4 Fan, J. S.; Zhang, M. Neurosignals 2002, 11, 315-321. 

5 Bezprozvanny, I.; Maximov, A. Froci. Natl. Acad. Sci. 
USA 2001, 98, 787-789. 

6 Kachinsky, A. M.; Froehner, S. C; Milgram, S. L. J. Cell. 
Biol. 1999, 145, 391-402. 



8 



7 Faulkner, G.; Pallavicini, A.; Formentin, E.; Comelli, A.; 
Ievolella, C; Trevisan, S.; Bortoletto, G.; Scanna- 
pieco, P.; Salamon, M.; Mouly, V.; Valle, G.; Lan- 
franchi, G. J. Cell. Biol. 1999, 146, 465-475. 

8 Sheng, M.; Sala, C. Annu. Rev. Neurosci. 2001, 24, 1-29. 

9 Hung, A. Y.; Sheng, M. J Biol Chem 2002, 277, 5699- 
5702. 

10 Doyle, D. A.; Lee, A.; Lewis, J.; Kim, E.; Sheng, M.; 
MacKinnon, R. Cell 1996, 85, 1067-1076. 

11 Au, Y.; Atkinson, A. R.; Guerrini, R.; Kelly, G.; 
Joseph, C.; Martin, S. R.; Muskett, F. W.; Pallavicini, A.; 
Faulkner, G.; Pastore, A. Structure 2004, 12, 611-622. 

12 Hillier, B. J.; Christopherson, K. S.; Prehoda, K. E.; 
Bredt, D. S.; Lim, W. A. Science 1999, 284, 812-815. 

13 Fujii, N.; Haresco, J. J.; Novak, K. A.; Kuntz, I. D.; 
Guy, R. K. Molecular Biology of the Cell 2002, 13, 360A- 
360A. 

14 Fujii, N.; Haresco, J. J.; Novak, K. A. P.; Stokoe, D.; 
Kuntz, I. D.; Guy, R. K. J. Am. Chem. Soc. 2003, 125, 
12074-12075. 

15 Shan, J. F.; Shi, D. L.; Wang, J. M.; Zheng, J. Biochem- 
istry 2005, 44, 15495-15503. 

16 Saro, D.; Li, T.; Klosi, E.; Udugamasooriya, G.; 
Spaller, M. R. Biopolymers 2003, 71, 413-413. 

17 Wiedemann, U.; Boisguerin, P.; Leben, R.; Leitner, D.; 
Krausc, G.; Moelling, K.; Volkmer-Engert, R.; Oschki- 
nat, H. J. Mol. Biol. 2004, 343, 703-718. 

18 Udugamasooriya, G.; Saro, D.; Spaller, M. R. Organic 
Letters 2005, 7, 1203-1206. 

19 Stiffler, M. A.; Grantcharova, V. P.; Sevecka, M.; 
MacBeath, G. J. Am. Chem. Soc. 2006, 128, 5913-5922. 

20 De Los Rios, P.; Cecconi, F.; Pretre, A.; Dietler, G.; 
Michielin, O.; Piazza, F.; Juanico, B. Biophys. J. 2005, 
89, 14-21. 

21 Tama, F.; Sanejouand, Y.-H. Prot. Eng. 2001, 14, 1-6. 

22 Delarue, M.; Sanejouand, Y.-H. J. Mol. Biol. 2002, 320, 



1011-1024. 

23 Gerstein, M.; Krebs, W. Nucl. Acid. Res. 1998, 26, 4280- 
4290. 

24 Madsen, K. L.; Beuming, T.; Niv, M. Y.; Chang, C. W.; 
Dev, K. K.; Weinstein, H.; Gether, U. J Biol Chem 2005, 
280, 20539-20548. 

25 Basdevant, N.; Weinstein, H.; Ceruso, M. J Am Chem 
Soc 2006, 128, 12766-12777. 

26 Gianni, S.; Engstrom, A.; Larsson, M.; Calosci, N.; 
Malatesta, F.; Eklund, L.; Ngang, C. C; Travaglini- 
Allocatelli, C; Jemth, P. J. Biol. Chem. 2005, 280, 34805- 
34812. 

27 Songyang, Z.; Fanning, A. S.; Fu, C; Xu, J.; Marfa- 
tia, S. M.; Chishti, A. H.; Crompton, A.; Chan, A. C; 
Anderson, J. M.; Cantley, L. C. Science 1997, 275, 73-77. 

28 Clementi, C; Nymeyer, H.; Onuchic, J. J. Mol. Biol. 
2000, 298, 937-953(17). 

29 Levy, Y.; Wolynes, P. G.; Onuchic, J. N. Proc Natl Am 
Soc USA 2004, 101, 511-516. 

30 Cecconi, F.; Micheletti, C; Carloni, P.; A., M. Proteins: 
Structure, Function, and Genetics 2001, 43,. 

31 Gianni, S.; Calosci, N.; Aelen, J. M.; Vuister, G. W.; 
Brunori, M.; Travaglini-Allocatelli, C. Protein Eng Des 
Sel 2005, 18, 389-395. 

32 Morriss, G. P.; Dettmann, C. P. CHAOS 1998, 8, 321-336. 

33 Ferrcnbcrg, A. M.; Swendsen, R. H. Phys. Rev. Lett. 1988, 
61, 2635-2638. 

34 Valleau, J. P.; Card, D. N. Journal of Chemical Physics 
1972, 57, 5457-5462. 

35 Gardiner, J. Handbook of Stochastic Methods; Springer- 
Verlag: 2002. 

36 Mel'nikov, V. I.; Meshkov, S. V. The Journal of Chemical 
Physics 1986, 85, 1018-1027. 

37 Klimov, D. K.; Thirumalai, D. Phys. Rev. Lett. 1997, 79, 
317-320. 



